Mach cones in the quark-gluon plasma: Viscosity, speed of sound, and effects of finite 

source structure 
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I use the space-time distribution of energy and momentum deposited by a fast parton traversing 
a perturbative quark-gluon plasma as a source term for the linearized hydrodynamical equations 
of the medium. A method of solution for the medium response is presented in detail. Numerical 
results are given for different values of the shear viscosity to entropy density ratio, 77/s, and speed 
of sound, Cs. Furthermore, I investigate the relevance of finite source structure by expanding the 
source term up to first order in gradients of a 5 function centered at the fast parton and comparing 
the resulting dynamics to that obtained with the full source. It is found that, for the source term 
used here, the medium response is sensitive to the finite source structure up to distances of several 
fm from the source parton. 
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I. INTRODUCTION 

A relatively new and exciting problem in quark-gluon 
plasma (QGP) physics is to determine the response of 
the medium to the passage of a fast parton. Fast partons 
are created by hard transverse scattering in the early 
moments of a heavy-ion collision and have long been 
considered a useful probe in understanding the QGP. 
The primary emphasis has focused on the process of jet 
quenching in which fast partons lose energy and mo- 
mentum by interacting with the surrounding medium 
(see, for instance, Rcfs. [H, i, i, i, S @, 0, S)- Re- 
cently, the question of how the energy and momentum 
deposited by the fast parton affects the bulk behavior of 
an evolvin g Q GP has gained attention (see, e.g., Refs. 
i, E, EH, III [II, H, [il, [H, [13, El)- Interest in un- 
derstanding the medium's response to the passage of a 
fast parton has been spurred on by experimental mea- 
surements at the Relativistic Heavy-Ion Collider (RHIC) 
[Tg . I20I l2l| of hadron correlation functions which suggest 
the fast parton may produce a propagating Mach cone in 



Distance scales in weakly coupled QGP 



Debye 



Boltzmann 1 Hydro- 



(g T)-^ 



(g'T)' 



FIG. 1: Some of the distance scales relevant to the hydrody- 
namic response of a weakly coupled QGP to fast partons. 
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the medium. 

There is strong evidence [22, ^M, that the matter pro- 
duced at RHIC obeys the hydrodynamic assumption of 
local thermal equilibrium. For this reason, the common 
theoretical approach to examining the QGP's response 
to a fast parton has been to treat it as a source of energy 
and momentum coupled to the hydrodynamic equations 
of the medium. This makes sense provided the medium 
maintains local thermal equilibrium following the passage 
of a fast parton. Assuming the medium does respond 
hydrodynamically to a fast parton then raises the ques- 
tion of what the distribution of energy and momentum 
deposited is. It has been observed [3, [ll|, [l3 that the 
medium's response to fast partons is sensitive to the spe- 
cific form of energy and momentum deposition, creating 
the need for a hydrodynamic source term derived from 
first principles. 

It is instructive to consider the mechanism of energy 
and momentum deposition and the different scales in- 
volved. In a quantum chromodynamic (QCD) plasma 
fast partons interact with the medium at a distance scale 
of the order of the inverse Debye mass, {mn)^^ . This 
interaction creates a disturbance, which in turn inter- 
acts with the surrounding medium, creating a new dis- 
turbance at some larger distance scale. The new distur- 
bance again interacts with the surrounding medium, and 
eventually the initial disturbance propagates outward to 
some arbitrarily large distance scale. At distances much 
greater than the mean free path. Ay, the medium's re- 
sponse to the initial disturbance can be accurately de- 
scribed by hydrodynamics. Thus an effective QCD hy- 
drodynamic source term should include the medium's re- 
sponse up to distances of order Ay, at which point the 
system evolves hydrodynamically. Whether the plasma 
is weakly or strongly coupled, the initial energy and mo- 
mentum deposition occurs at a distance scale of the order 
of the inverse Debye mass, although the specific value of 
the Debye mass depends on the strength of the coupling. 
However, in a strongly coupled plasma, the concept of a 
mean free path loses meaning, and instead the de Broglie 



wavelength sets the minimum scale at which the hydrody- 
namical description is valid. In general, the application 
of hydrodynamics is valid on shorter distance scales for 
more strongly coupled mediums. 

Recently, Neufeld [25 presented a derivation of the 
hydrodynamic source term expected from a fast par- 
ton moving through a pcrturbativc QGP, both with and 
without including the effect of color screening, by includ- 
ing the medium response at a distance scale of the or- 
der (771 £))^^. Using the unscreened, relativistic form of 
this source term coupled to the linearized hydrodynam- 
ical equations of the medium the authors of Ref. |18| 
showed that the medium response includes a propagat- 
ing sound wave with the shape of a Mach cone and a 
diffusive wake. In this work, I will use a slightly modified 
form of the relativistic limit of this source term in the 
linearized hydrodynamical equations of the medium. A 
detailed solution of the equations of motion will be pre- 
sented along with the resulting dynamics for a range of 
values of the shear viscosity to entropy density ratio, yy/s, 
and speed of sound, Cg. I will also expand the source term 
up to first order in a series of gradients of a, S function 
and compare to the full result, in an effort to understand 
the relevance of finite source structure. It will be shown 
that the medium response is sensitive to the finite source 
structure up to distances of several fm from the source 
parton for the source term used here. 

In a weakly coupled QCD plasma at high temperature 
T, the inverse Debye mass is of order {gT)~^, where g 
is the running coupling, whereas the transport mean free 
path is of order {g^T)~^ [2J]. The medium's response 
to disturbances at distance scales between {gT)~-^ and 
(g^T)^^ is accurately described by the Boltzmann equa- 
tion (see Fig. [T|). The hydrodynamic source term ex- 
amined here includes the medium's response at a dis- 
tance scale of the order of the inverse Debye mass, at 
which point the medium's hydrodynamic response is in- 
voked. This is a simplification of the true QCD evolu- 
tion, where the medium's evolution between (gT)~^ and 
(f/^T)^^ should be described by the Boltzmann equation, 
after which the hydrodynamic response can be invoked. 
However, in the QGP produced at RHIC, it is likely that 
the mean free path is comparable in size to the inverse 
Debye mass (compare, for instance, Refs. [2a] and [27|). 
Thus, from a phenomenological point of view, physics at 
a distance scale of the order (tti^i)"^ may be relevant to 
the specific structure of the medium's hydrodynamic re- 
sponse to fast partons. However, this also suggests the 
QGP produced at RHIC may be strongly coupled, limit- 
ing the application of perturbation theory. 

The paper is structured as follows. In SeclHJ I con- 
sider the general form of linearized hydrodynamics with 
a source term. I then introduce the specific source term 
studied here. In Sec lIIII the source term is expanded in 
terms of gradients of a J function up to first order. Both 
the full source term and the d function expanded form 
arc then Fourier transformed into momentum space. In 
Sec lIVI the hydrodynamic equations are solved in terms 



of a one-dimensional numerical integration. In Sec(V| re- 
sults are given for both the full source and the 5 function 
expanded version, as well as different values of the vis- 
cosity and speed of sound. In what follows I choose units 
such that h = c = kt = I. 



II. LINEARIZED HYDRODYNAMICS WITH A 
SOURCE TERM 

The first-order hydrodynamical equations for a 
medium with nonzero shear viscosity rj in the presence 
of a source term J'' are given by 



gJ.^.. ^ ju^ 



(1) 



where T^'^ is the energy- momentum tensor of the system. 
If one assumes that the energy and momentum density 
deposited by the fast parton is small compared to the 
equilibrium energy density of the medium, the hydrody- 
namical equations [Eq. ([1])] can be linearized. Defining 
TA"- ^ T^^" + STf"", where STf"" is the perturbation of 
the energy-momentum tensor resulting from the source 
in an otherwise static medium, one has 

a^jT'"' = J^ (2) 

where d^.T^'' = and ^T'"' is given by [I 



Sr^ = 6,jclSe - -F, d'g^ + d^ g' - -Jy V • g 



(3) 



In 



Eqs. © Cs denotes the speed of sound, F 



47; 



^T7 — , — T — T^ is the sound attenuation length, and en 
and po are the unperturbed energy density and pressure, 
respectively. 

By introducing the general rule for Fourier transforms 

the equations given by Eq. ([2|) arc written in momentum 
space as 



J° = —iujSe + ik • g 



(5) 



J = -tu;g + ikcide+^rJkh+^{k-s)Y (6) 

Solving for k • g in Eq. ([6|) allows for Se to be determined 
fromEq. ©: 



k • g = kgL 



k • J - ik'^c^Se 

—iuj -\- Fj/c^ 



(7) 



and hence 






where the source and perturbed momentum density vec- 
tors are divided into transverse and longitudinal parts: 
g = kgL + gT and J = kJ^ + Jt, with k denoting the 
unit vector in the direction of k. Similarly, one has from 
Eqs. © and dS]) 



kgL = —iJ + (jjSe 



yielding 



gL(k,w) =k5L 



iLukJL{k, uj) + ic2k,/"(k, LU) 
Lu^ — c^k^ + iTsUjk^ 



(9) 



(10) 



The transverse part of g can be obtained from Eq. ^. 
The calculation is simplified by considering that any part 
of g proportional to k is a part of gi. This leaves 



gT(k,w) = g-gL 



iJT(k, uj) 



(11) 



Equation (|lip is a diffusion equation and the quantity 
gT is interpreted as diffusive momentum density gener- 
ated by the fast parton. Equations ([8|) and pO]) describe 
damped sound waves propagating at speed c^: it follows 
that Se and g^ are interpreted as the energy and momen- 
tum density carried by sound generated by the fast par- 
ton. The importance of the explicit form of the source 
term can be readily seen. In a homogeneous medium 
symmetries ensure that the source vector, J, can be writ- 
ten generally as 



J(x,t) = up{r)+Vq{r) 



(12) 



where u is the velocity of the source particle and p{r) 
and q{r) are scalar functions of the (possibly Lorentz- 
contracted) magnitude r = •\/(x — ut)^. If we instead 
write Eq. (J12p in momentum space we have 



J(k, w) = u / d*x e*'=-"p(r) + ik / d^x e'^-^'qir) (13) 



ik-x , 



and 3t is found to be 

k(k- J) 



JT(k, w) = J - 



fc2 
Uk^ - k(k • u) 



(14) 



(fxe"'-''p{r) 



so that q{r) does not contribute to Jt- If one chooses 
the source [Eq. (fT2|) ] such that p{r) = then there is no 
excitation of the diffusive momentum density. It is clear 
that the hydrodynamics of the system are sensitive to 
the specific form of the source term. In particular, terms 
that are in the form of a gradient only generate sound. 

As mentioned previo usly , in this work I use the source 
term derived in Ref. [2^. There, the fast parton was 
treated as the source of an external color field interacting 
with a perturbative QGP through a Vlasov equation. For 
a gluonic medium at temperature T in the presence of a 
parton moving with velocity u = uz at position r — utz 



in the relativistic limit (7 = {1 — v?) ^1"^ ^ 1), the source 
is given by 



J^{x) = (J"(x),uj*'(a;)-Jv) 



where 



J°(p,z,f) = d{p,z,t)ju^ 1- 



3v{p,z,t) = {x— ut) d{p, z,t) 



^UZ- 



Z^ 72 + p2 



2:^.72 + p2 



(15) 

(16) 
(17) 



and 



d{p,z,t) 



87r(p2+ ^2^2)3/2 ■ 



(18) 



In the above expressions, (Qp)^ = 3 for a gluon and 4/3 
for a quark, p = (x^+y^)^' ^ is the radius transverse to the 
z axis, as ~ .<7^/47r is the strong coupling, tod = gT and 
z_ = (z — ut). In what follows any numerical coefficient 
suppressed by powers of 72 will be dropped. For instance, 
terms such as 72 -|- 1 will be taken as 7^. 

The vector part of the source, J = uJ° — Jv, is explic- 
itly 



asjQ^rmlu 
Btt 



M^ (a:,y,Z-72) 

(^2^2.^^2)3/2 (^2^2.,_^2)2 



7U 



(19) 



which, as one can verify, can be re-written in the form of 



Stt 



7U , ^_ u^ 



(zi72+p2)3/2 2(2^72 +p2) 

(20) 



III. DELTA FUNCTION EXPANSION AND 

MOMENTUM SPACE REPRESENTATION OF 

THE SOURCE 

At distances increasingly far from the fast parton, one 
expects that the source term will begin to look like a 6 
function. Since hydrodynamics is a long distance effec- 
tive theory, the hydrodynamic solutions (in the range of 
validity) should be dominated by the lowest order terms 
in an expansion of gradients of a (5 function centered at 
the location of the source parton (a detailed discussion 
of this is given in [l6|). Higher order terms in the ex- 
pansion, which are sensitive to the detailed structure of 
the source term, will become less important at larger dis- 
tances. In this section, I will expand the source term, as 
given by P^ and (PO)) . up to first order in gradients of 
a 6 function. Later, the hydrodynamic equations, ([H]), 
(fTU)) . (fTTj) . will be solved for both the full source term 



and the truncated scries. A comparison of the solutions 
will highlight at what distance scales the detailed struc- 
ture of the source term becomes negligible. It will prove 
convenient to Fourier transform the source into momen- 
tum space, which I will also do in this section, before 
attempting to solve the hydrodynamic equations. The 
effect of color screening, which is absent in (fTB|) and (PO)) . 
will be modelled by including a damping factor of the 
form e"'"", where m^^ is a typical screening scale. In a 
perturbative QGP the inverse screening scale is given by 
TTiD = gT, which appears as a coefhcient in front of the 
source term used here. However, in principle, at higher 
orders the screening scale may be different than what 
appears as the coefficient of the source term. It is thus 
instructive to keep m arbitrary; however, in solving the 
hydrodynamical equations in Scc lIVI I will set m = gT . 
Also, when necessary, a short distance cutoff will be used 
to regulate ultraviolet divergences. A common choice, 
which will also be used here, for the short distance cutoff 
in colhsional energy loss is Pmin = {'^\/ EpT)~^ , where Ep 
is the energy of the fast parton (sec, for instance, [3(J|). 
Consider (|16p . which can be expanded as 



J\p,z,t) = Co(5(x_) + Ci • V<5(x_; 



where I have used the shorthand notation 



5{iL_) = 5{x)5{y)5{z - ut). 



(21) 



(22) 



The coefficients, Cq and Ci, arc found by taking the ap- 
propriate moment of J''(/9, z, t). Introducing the damping 
factor, e~'"^, one has for Cq, 



Co= / d^xJ^{p,z,t)e-P 
(fx 



asiQ^pfml-^u^ ^_^^^ 



87r(p2+ 72^2)3/2 

aN2^2 „,2 / „ 



(23) 



7. "^0 



2/^r 



where {2y/EpT) ^ has been introduced as a short dis- 
tance cutoff, and Go is a representation of the incomplete 



Gamma function 



Go{z) 



00 g-t 

dt . 

t 



(24) 



The coefficient given by ([25)1 gives the total energy de- 
posited into the medium per unit time. 
Similarly, Ci can be obtained as 



Ci 



= / dh 



d^x{x,y,z_)J^{p,z,t)e-f"^ 

asjQlfml^ ^u^ z_{x,y,zJ) 
Stt 



C^,,(Qg)2^2^^2 



0,0,- 



(p2 -1-72^2 )2 



-p in 



(25) 



2 V ' '4m7^ 

The results from ((231) and ([25]), together with jUl), give 



Jl{p.z,t)^ 



asiQpYmf^u 



(26) 

where the subscript d is meant to indicate the expansion 
in gradients of a 5 function. Proceeding in an analogous 
manner yields for (P0|) 



Jr,(/?,z,t) = i^- uGo 



2^Epf 



(27) 



+- (w^v + u(u • V)) (5(x_) + 

8 m7 / 



Equations ^^ and ^7} provide the expansion of the full 
source, (|T6|) and ((20)l . up to first order in gradients of a 
d function. 

As previously mentioned, it is easiest to solve for the 
hydrodynamics in Fourier space. To do this, it is neces- 
sary to first transform the source terms into momentum 
space, following the general rule (JH). The details of the 
Fourier transforms of p6p and (j20p are give in Appendix 
A, and the result is 



J°(k,c.) = 
J(k,c.) = 



asiQp) m^u 



■{2tt)S{lu - uk:,) Go 



as{Q;fmlu^ 



{2TT)d{uJ - Uk,) 



m + fc-r 



i7r(u • k) 



2/EpT/ 4-f^l4T: 



(28) 



X uGo 



fc7 



iiT I \[W^ 



2/EpT/ 47 






,, u(u • k)m 



VM 



These equations should be compared to the Fourier transforms of (|26p and (j27p . which arc found by making the 



replacements S{z-) -^ (27r)5((jj — uk^) and V -^ ik: 



J°,{k,u;) = '^^^'^ " [2i,)5{u - uk,) f Go 




2,/E^ I 8m7 



(29) 



k + u(u-k)) 



r 



One can verify that (pS)) reduces to ([M)) by taking the 
kx ^ hmit in the coefficients of 1 and k. (P5|) and (P5|) 
wiU be used in the next section to solve for the hydrody- 
namic variables of the medium. 

As discussed at the beginning of the section, at dis- 
tances increasingly far from the fast parton the source 
term is dominated by the lowest order terms in an expan- 
sion of gradients of a J function. It's clear from inspec- 
tion that the detailed structure of the full source term, 
given by (j28p . becomes important at a momentum scale 
/c ~ m. This could have been anticipated, since the full 
source term is calculated up to distances of the order of 
the screening length. A quantitative comparison of the 



effects of psp and ([29]) requires solving the hydrodynamic 
equations. This will be done in the next section. 



IV. SOLVING THE EQUATIONS 



The result given by ([28|) is combined with equa- 
tions dHD, Cni), and (HID to yield 5e(k,a;), gL(k,a;) and 
gT(k, w). These are then transformed back to position 
space using the relation (j4]). In Se(k,uj) and gL(k, w), one 
can find J^ by taking k • J, which is written conveniently 
as 



J 



Jl = "^L i27r)S{uj - uk. 



2fc 



(u-k)jo as{QlYm'j,u 



(u • k) Go 



mo 



+ '•'":. " {2t:)5{oj - uk.Y-^^^^ ( ^|k? ' -2 



2k 



47 



i7r(u • k) 



2/EpT / 47^^ 



■^f, + mj, — mD 



47 



k'^ + ?71q — TTlD 



(30) 



where I am now taking m ~ mu- After integrating out 5{uj ~ uk-,) and using (jA-Sp . the expression for (5e(x, t) is given 
by 



(5e(x,i) = 



87r2c2 
mo + fcr 



Go 



2/EpT / 47y^|~r 






(2iuk,-T,k^) - 



47 






(31) 



^D - mu 



where A^ = c^/i^^ ^ c^), a = Tsu{X^ /c1)kz{k^ + k1) and I am again working in plane polar coordinates. The integral 
over kz can be performed using contour integration. Poles are located at kz = ±(fc|,A^ =F ijcrl)^/^, where \a\ is itself 
a function of k^. When evaluating the residues at these poles I make the approximation a{kz) ~ cr(±fcyA). This 
approximation is valid at momentum scales for which the sound attenuation is small (fcy <C Cj/Fs), which should be 
reasonable in the hydrodynamic limit. Both poles are located in the lower complex plane so that the integration only 
contributes for z < ut, i. e., behind the source parton. Performing the integration yields 



..(x..)^°-"^g!'":°"'^'R. 




dh 



iJoipkT)e 



ikTV ^^ ~'^f^Taiz—ut) 



y/X^ — ikxa 



TTU 

47 






If) - m.D 



iirukT^/ X'^ — ikTOi 



I 2iukT\/X'^ — ik^a — Tskrpll + A — ikxd) ) 



(32) 



where a = V suX"^ / {c^^{\^ + 1)). The final integration over fc-r is performed numerically. The analogous expression 
resulting from ([^5]) is given by 



(5eD(x, t) 



Q,(Qg)^mg,7.^A^ 



Go 



I mo 



Re 



dh 



iJo(pfcT)e*''^^^'"'''^"^^""*^ /TTW^fc, 



VA^ — ikTO: 



,2l2 



877710 



2/^ 



iTrukT\/X^ — ikxa 
477710 



f 2iukT\/X^ — ikxa — rsfcy(l + A — ikTo) 



(33) 



The same approach is applied to g^ (x, t) . The contour integration proceeds in the same manner as in (j3ip with 
the exception that one term has additional poles at k^ = Hkx- The additional pole at kz = ikx allows for some 
contribution in the region in front of the source parton. The result is 



for z < ut : 
gL(x,i) = '- — Re 



47rc? 



dkT .,;. -iJiipkx), -iJiipkr), V A^ - ikraJoipkr) x 



Go 



mo + kT\ f 7i^(A^ — ikro) 



c + 



iTTukT y/X^ — ikra 



2y/E^ J \{l + X^-ikTa) ' ""J ' A-/ y^j, ^ , , .-^ 

y'dfcTe-'^-l^^"*! (--MpkT),-^MpkT),MpkT)) 



{u' + cl) + 



^■mD{mr, - a/Zcj + 777,p) 
fc|,(l + A2 - ikra) 



a,{Q;)^mlu^ 



X kr Go , + — — 1 

V v^v/^; 47 1^ 

for z > ut : 

asiQtfmlu^ 



rriD 



Vk'. 



y + 777,^ 



gL(x,i) =- 
X fcr Go 



mjj + kT \ TTnio 



dkr e-'=-l^-"*l (-Mpkr), ^Mpkr), Mpkr) 
\P P 



1- 



771D 



2/Epr/ 47 \ yfc2 



which must be done numerically. The analogous expression resulting from p9p is given by 



(34) 



for z < ut 
gL^(x,t) = 

I Go 



a,{Q;)^mlu^X' 



Re 



47rc2 
mo \ / 77^ (A^ — ikxa) 



jtj, ^ikTV^'-'-ikraiz-ut) /^ y , 

dkT ,,„ -iJiipkr), -iJiipkr), V X^ - ikraJoipkr) x 

VX^-ikra \P P 



2y%fJ \il + X^-ikTa) 



c: + 



47 777D 



(7.^+C^) + 



+ 

for z > ut 
gL^(x,t) = - 



2(1 + A2 - ikra) 
mo \ . 



Trfc^ 



"^^"^^^ "'"^ /dfc^e-'^-l— *l f-^Ji(pfc^),-^Ji(pfc^), Jo(pm) (fcTG, . /_ , , ^ 

St: J \ p p J y y2^EpT J 8"/m-D 

as{Q;ymlu^ 



Stt 



dfcTe-'=-l^-"*l r-Ji(/7fcT),-Jl(pfcT), Jo(pfcT)) (fcTGo 



777D 



Trfc^ 



2y/E^ 8777115 



(35) 



Combining (fTil) with (fTTI) and following the same approach used above allows for the determination of gT(x, t). 
The denominator of (|lip has a simpler structure than in (jSllOp which allows the contour integration to be performed 



exactly. Defining il = 4u/3rs it is found that 



gT(x,f) = f- / dkr- 



4tt 



(4-(ZiT)2)Wl + ^ 



4feJ 



16u7r 
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J \ P P 



(36) 



^ 2/EpT y ' 47 

where =F refers to the sign of [z — ut) and 



fviji \/ i\"-p ~r ^^T) 



-4'-y-f 



(37) 



The analogous result from (|29p is given by 
gToi^^t) = 7 / dkr- 



Att 



ikl-{ATy)Ji + ^ 



Go 
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2/EpT/ 877710 / V P P 



± 



asiQ^rml 
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NUMERICAL RESULTS AND DISCUSSION 



Having obtained expressions for the hydrodynamic 
quantities Se{x,t), gL(x, f) and gT(x, f) I now consider 
the results of numerical integration. All calculations are 
performed for a gluon moving along the positive z axis 
at position ut and speed u = 0.99955c (7 w 33). The 
strong coupling, as, is chosen to be I/tt, the temperature 
is taken to be T = 350 MeV, and Ep = 16 GeV. As men- 
tioned in the introduction, I will compare the solutions 
resulting from the full source term and the S function ex- 
panded source term, given by ((28)) and ((29l) , respectively. 
I will also make a comparison of the results for a range 
of values of the shear viscosity to entropy density ratio, 
rj/s, and speed of sound, c^. 

The solutions resulting from the the full source term 
and the S function source are shown in Figures [2] and [31 
In both Figures the results are plotted for 77/s = l/47r 
and Cs = c/VS. Figure [2{a) shows a contour plot of 
(5e(x, t) for each source. In both cases a well define Mach 
cone is visible in the trailing medium. Figure [2Kb) shows 
the result for the magnitude of the momentum density, 
|g| = ISLiy^jt) + gT(x, t)|. One now sees both a sound 



contribution from gi (x, t) , which excites a Mach cone, 
and a diffusive contribution from gy(x, t), which is ex- 
cited in the region directly behind the source gluon. The 
diffusive momentum density produces flow almost exclu- 
sively in the direction of the source parton's velocity, 
while the Mach cone generates flow outward and perpen- 
dicular to it's boundary, as indicated by the arrows. One 
can see from Figure [5] that the full source term and the 6 
function expanded source term provide qualitatively sim- 
ilar results, particularly in the region far from the source 
parton. In the region near z_ = 0, the full source solu- 
tion has a noticeably larger transverse extent than the 
corresponding 6 source result. A more quantitative com- 
parison can be made by examining Figurc[3|where 5e(x, t) 
is plotted as a function of z- for fixed p. Here one sees 
that the two results converge at a distance of about 5-6 
fm behind the source parton for the chosen p values. 

Results arc next presented for three different values of 
the shear viscosity to entropy density ratio, rj/s. The first 
value chosen for rj/s is l/47r « 0.08, which has been pro- 
posed [2^ as a universal lower bound for all relativistic 
quantum field theories and is calculated in the strongly 
coupled limit. The other two values for rj/s are multi- 
ples of the previous value, 3/47r and 6/47r, and are more 
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FIG. 2: (Color online) Plots of (a): the perturbed energy density, and (b): the perturbed momentum density, for both the 
full source term and the S function expanded source term. Here, rj/s = l/Air and Cs = 0.57. The direction of the momentum 
density is indicated by the arrows. 
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FIG. 3; (Color online) A comparison of the perturbed energy density generated by the full source term and the & function 
expanded source term for two different values of p. The two results converge at a distance of about 5-6 fm behind the source 
parton for the chosen values of p. 



consistent with the application of perturbation theory, 
which is the method used to calculate the source used 
in this paper. For example, Arnold el al. [29[ found for 
the leading order result rj/s = 0.48 for a gluonic plasma 
with as = 0.3. More recently, Xu and Grcincr found 



jy/s = 0.13 for a gluonic plasma with the same value of 
as by going beyond leading order in the diluteness of 
the medium [3l|. A small value of the shear viscosity, 
which is required by the RHIC data |32l | , is not necessar- 
ily incompatible with perturbation theory, especially if 
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FIG. 4: (Color online) Plots of (a): the perturbed energy density, and (b): the perturbed niomentum density for different 
values of the shear viscosity to entropy density ratio, r\l s. The black lines in (a) are drawn where one would expect to find 
the boundary of a Mach cone in the absence of dissipative effects. Plots scaled by the radius, p, which factor in the conical 
broadening of the cone, are shown in Figure [S] 
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FIG. 5: (Color online) Plots of (a): the total perturbed energy density, and (b): the total perturbed momentum density, 
contained at a given radius in the p — z- plane for different values of 77/s. As one can see in (b) the total perturbed momentum 
density carried by the sonic Mach cone exceeds that contained in the diffusive wake. 
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FIG. 6: (Color online) Plots of (a): the perturbed energy density, and (b): the perturbed momentum density for different 
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in Figure [T] 



2jrp X & [GeV/fai''] for cs = 0.57 



cs = 0.45 



cs = 0.3 




[fm] 



2jrp X Igl [GeV/fin''] for cs = 0.57 




5j2_ [fm] 




cs = 0.3 



[fm] 




5s_ [fm] 



p[fm] 



FIG. 7: (Color online) Plots of (a): the total perturbed energy density, and (b): the total perturbed momentum density, 
contained at a given radius in the p — z- plane for different values of Cs . The magnitude of the Mach cone, when integrated, 
is similar for the different plots. 



11 



the viscosity is lowered by anomalous contributions [33| . 

The results for (5e(x, t) and |g| = |gL(x, t) +gT(x, t)| 
for all three viscosities are shown in Figures S] and [51 
where Cg — c/v3- The black lines in the contour plots 
of Figure [H^a) are drawn along the slope x = ±\{z — ut) 
which is where one would expect to find the boundary of a 
Mach cone in the absence of dissipative effects. In Figure 
[5l the total energy density (figure (a)), and magnitude 
of momentum density (figure (b)), contained at a given 
radius are shown. It is clear from the plots that the Mach 
cone broadens and weakens as the viscosity is increased. 

Finally, results are shown for three different values of 
the speed of sound, c^. The first value is c^ — 0.57, which 
is the the limiting value for a conformal ideal relativistic 
gas, while the other two values arc c,, = 0.45,0.3. It's 
likely that the QGP produced at RHIC experiences a 
speed of sound close to all three of these values during 
its evolution [3J|. The results are shown in Figures [5] and 
[71 where I have chosen 77/s = I/A-k. One should note that 
the diffusive contribution, gy (x, t) , is independent of the 
speed of sound. 

It is interesting to consider how the results presented 
here compare to experimental data. In the di-hadron cor- 
relation functions measured at RHIC there is a double 
peak structure in the back-jet (source parton) distribu- 
tion which has been interpreted by some as the result of 
Mach cone generated flow. In the spectrum presented 
here one indeed finds Mach cone generated flow but also 
finds a substantial diffusive flow, which seems to be miss- 
ing from the RHIC data. In order to make a comparison, 
however, one must consider that the matter created in 
heavy ion collisions at RHIC rapidly expands in contrast 
to the static background assumed here. The diffusive 
momentum is deposited locally and is thus probably dif- 
ficult to observe in an expanding medium. On the other 
hand, the Mach cone propagates at the speed of sound, 
which is of the same order of magnitude as the expansion 
velocities in the matter produced at RHIC, and is likely 
more readily observed experimentally. 

Care should be taken when examining the azimuthal 
particle spectrum generated by a fast parton using 
an isochronous Cooper-Frye freeze-out from a static 
medium, such as in the work done by Betz et al. [35J . 
In their paper, the authors compared the perturba- 
tive QCD based source term studied here with one de- 
rived in the strongly coupled Ads/CFT correspondence 
[16| . Their conclusion was that the anomalous azimuthal 
hadron correlations observed at RHIC are likely the re- 
sult of flow generated by the non-equilibrium Neck zone 
in the Ads/CFT case, a contribution which does not obey 
Mach's law [30] • In the isochronous Cooper-Frye freeze- 
out scenario the entire volume of matter is assumed to 
hadronize at the same time, independent of physical pro- 
cesses. This is in contrast to the freeze-out in a heavy ion 
collision, which occurs as the result of an expanding and 
cooling medium. The effect of an isochronous freeze-out 
scenario is that any cylindrically symmetric, or conical, 
contributions tend to be washed out (see the discussion 



in [37[). Any rigorous comparison to experimental re- 
sults will require incorporating a realistic source term in 
an expanding medium. 

In summary, I have here presented in detail a method 
of solution for the linearized hydrodynamical equations 
of a QGP coupled to the source term generated by a 
fast parton. The solution has been examined for differ- 
ent values of the shear viscosity to entropy density ratio 
and speed of sound. Additionally, the relevance of finite 
source structure has been investigated by performing an 
expansion in gradients of a 5 function centered at the lo- 
cation of the source parton. Comparison of the medium 
response generated by the full source with that generated 
by the 5 expanded one shows that the result is sensitive 
to the finite structure up to distances of several fm from 
the fast parton for the source examined here. 
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Appendix A: Obtaining the Fourier Representation 
of ([16]) and ([20l) 

In Sec lIIII the explicit determination of the Fourier 
transform of the full source was put off to this appendix. 
Including the damping factor, e~f™. one has for the in- 
verse Fourier transforms of ([TBI) ^^^^ (EOj) 
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After an integration by parts, the second term in (|A-1|) 
takes the form 



d^x 



y2 ^ik-x — pm 

2^,2 ! r.'2\ 



2{zlr 



(ik + m(cos(?!), sini/), 0)) (A-3) 



where I am working in plane polar coordinates, p and 0, 
such that X = p cos and y = p sin (j). 

It is clear there are three distinct integral forms which 
need to be evaluated. After (trivially) integrating out the 
t dependence to bring down a factor of 27r(5(a; — ufc^), the 
three distinct integral forms are 
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The exponential dependence upon the variable (/) is in 
the term —ip^k^ cos 4> + ky siiKJ)). Re- writing this term 
as —ipkx cos [(f) — a], where k^ = kxcosa and ky = 
/cysina, the (fi integration can be done using the rela- 
tions 
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where Ji{x) is the Bessel function of the first kind of 
order i. The final result for A2 is obtained by using the 
relation 
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To determine A3 it is necessary to evaluate 
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The above form of (|A-11[) is difficult to evaluate analyt- 
ically. However, it can be made more manageable by 
including the screening factor in the z, rather than p, 
integration. In this case, one has in the large 7 limit 
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where Zmin = i'^lx/EpT) ^ has been used as a short 
distance cutoff in the z integration. This gives for A3 



which gives, in the large 7 limit, 
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The first component of Ai, denoted by Ai^, requires 
evaluating 
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The second two components of Ai, denoted by Ai^, can 
be determined after evaluating 



dz / dp 

00 Jo 



e-'^'''-P'^Ji{pkT)p 



(p2 -)- 72^2) 



1- 



m7+ Ifc^l 



krj \ ^fc2 + {kj, + m2)72 + 2m7|fc2| 



(A-9) 



Again, working in the large 7 limit, this gives for Ai 
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It is now possible to write down the final result for 
(|A-2[) and (jA-ip . Remembering to include the factor of 
2t:S{lu — ukz) from the t integration, one has 
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(A-10) which is the result quoted in ([28 
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